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ABSTRACT 

We present a detailed analysis of the properties of tidally stripped material from dis- 
rupting substructure haloes or subhaloes in a sample of high resolution cosmological 
' r~| ' A^-body host haloes ranging from galaxy- to cluster-mass scales. We focus on devising 

pi |. methods to recover the infall mass and infall eccentricity of subhaloes from the prop- 

erties of their tidally stripped material (i.e. tidal streams). Our analysis reveals that 
there is a relation between the scatter of stream particles about the best-fit debris 
plane and the infall mass of the progenitor subhalo. This allows us to reconstruct the 
' infall mass from the spread of its tidal debris in space. We also find that the spread in 

radial velocities of the debris material (as measured by an observer located at the cen- 
tre of the host) correlates with the infall eccentricity of the subhalo, which allows us to 
reconstruct its orbital parameters. We devise an automated method to identify leading 
and trailing arms that can, in principle at least, be applied to observations of stellar 
streams from satellite galaxies. This method is based on the energy distribution of ma- 
terial in the tidal stream. Using this method, we show that the mass associated with 
. leading and trailing arms differ. While our analysis indicates that tidal streams can 

be used to recover certain properties of their progenitor subhaloes (and consequently 
satellites), we do not find strong correlations between host halo properties and stream 
properties. This likely reflects the complicated relationship between the stream and 
the host, which in a cosmological context is characterised by a complex mass accretion 
history, an asymmetric mass distribution and the abundance of substructure. Finally, 
we confirm that the so-called "backsplash" subhalo population is present not only in 
k>( \ galaxy cluster haloes but also in galaxy haloes. The orbits of backsplash subhaloes 

. brought them inside the virial radius of their host at some earlier time, but they now 

' reside in its outskirts at the present-day, beyond the virial radius. Both backsplash and 

bound subhaloes experience similar mass loss, but the contribution of the backsplash 
subhaloes to the overall tidal debris field is negligible. 

Key words: methods: n-body simulations - galaxies: haloes - galaxies: evolution - 
cosmology: theory - dark matter 
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1 INTRODUCTION 

The currently favoured model for cosmological structure for- 
mation is the ACDM model. This model assumes that we 
live in a spatially flat Universe whose matter content is dom- 
inated by non-baryonic Cold Dark Matter (CDM) and whose 
present-day expansion rate is accelerating, driven by Dark 
Energy (A). In the context of this model, the formation 
and evolution of cosmic structure proceeds in a "bottom- 
up" or hierarchical manner. Low mass gravitationally bound 
structures merge to form progressively more massive struc- 

* E-mail: kwarnick@aip.de 



tures, and it is through such a mergi ng hierarchy that 
galaxies, groups an d clusters form (e.g. iDavis et all Il985l : 
ISpringel et alTl2006l l. 

The tidal disruption of satellite galaxies around our 
Galaxy and others is characteristic of the hierarchical merg- 
ing scenario. As a satellite galaxy orbits within the gravi- 
tational potential of its more massive host, it is subject to 
a tidal field that may vary in both space and time. The 
gravitational force acting on the satellite strips a stream of 
tidal debris from it, and in some cases tidal forces may be 
sufficient to lead to the complete disruption of the system. 
This tidal debris tends to form two distinct arms - a leading 
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arm ahead of the sateUite and a trailing arm foUowing the 
sateUite. 

The process by which a sateUite galaxy suffers mass loss 
and the subsequent formation of a stream of tidal debris 
(hereafter referred to as a Udal stream) has been studied ex- 
tensively using numerical sim ulations (e.g. |johnstonlll998l : 
llbata et al]|200ll : iHelmil bool ) . A number of these studies 
have followed the orbital evolution of an individual satel- 
lite galaxy with a known mass profile, realised with many 
particles, in a static analytic and non-cosmo l ogicar host 



potential f e.g. I Johnston|[l998l : llbata^ et al.ll200ll : lHelmi 



2004 



and many others) . This approach allows for well-defined test 
cases to be investigated and for a systematic survey of orbital 
and st r uctur al parameters to be carried out. For example, 
iHelmil l|2004h showed that the width of the tidal debris in 
the plane perpendicular to its orbit increases when the host 
halo potential is fiattened rather than spherical. 

In this regard, it is important to note that stream 
properties are sensitive to both the internal properties of 
its progenitor satellite and the details of the orbit that it 
follows. For example, we expect the length and width of 
a stream to tightly correlate with the age of the stream 
and the mass/internal vel o city dispersion o f its progenitor 
(e.g. [johnston et al. 1 1 19961: |Johnstonlll998l : Ijohnston et al.l 
I2OOII : IPefiarrubia et al"]|2006h . Understanding what role the 
progenitor satellite and its orbit plays in shaping stream 
properties is therefore vital if we wish to use streams to 
address more fundamental problems, such as the nature of 
the dark matter. 

While previous studies have provided important in- 
sights into the formation and evolution of tidal streams 
and their dependence on the properties of both satellites 
and host dark matter haloes, none have considered tidal 
streams in cosmological host haloe^. Cosmological haloes 
have complex mass accretion histories and mass distribu- 
tions that can at best be crudely approximated by fitting 
formulae, and the interplay between a halo's asphericity, its 
dumpiness, and how its mass grows in time may be non- 
trivial. How this complexity may affect the results of pre- 
vious studies, which were based on analytic potentials and 
approximations to halo growth, is difficult to say. For ex- 
ample, interactions between the stream and substructures 
lead to the dynamical heating, causing an increase in inter- 
nal velo city dispersion and a subsequent broadening of the 
strea m (|Moore et al ] |l999l : llbata' et al ] |2002l : I Johnston et al.l 
|2002| ). Therefore, while studying tidal streams in cosmolog- 
ical haloes is a challenging problem, such an approach is 
necessary if our results are to refiect the complexity of na- 
ture (assuming the validity of the ACDM model, of course!). 

This paper is the first of in a series that will examine 
the formation and evolution of tidal streams in cosmological 
dark matter haloes, and determine how the properties of 



^ By "non-cosmological" we mean that the systems were evolved 
in isolation, in the absence of the background expansion of the 
Universe, without infall and accretion, and without a population 

of substructures. 

^ We note the work of IPefiarrubia et al.l 1I2OO6I ). which employs 
a semi-analytical approach to study tidal streams in an evolving 
host potential whose time variation is based on results of cosmo- 
logical simulations. 



these tidal streams depend on the properties of the satellite 
galaxies from which they originate and on the host dark 
matter haloes in which they orbit. We use high resolution 
cosmological simulations of individual galaxy- and cluster- 
mass dark matter haloes - the host haloes - that allow us to 
track in detail the evolution of a large population of their 
substructures or satellite galaxies over multiple orbit^. We 
seek to answer the questions we consider to be key to the 
study of tidal streams; 

• Can observations of tidal streams be used to infer the 
properties of their parent satellite? 

• Can tidal streams reveal properties of the underlying 
dark matter distribution of the host halo? 

The present paper primarily focuses on the first ques- 
tion and can be summarised as follows. In Section [2] we de- 
scribe the fully self-consistent cosmological simulations of 
structure formation in a ACDM universe that we have used 
in our study. The suite of host haloes (consisting of eight 
cluster and one Milky Way sized dark matter halo) are pre- 
sented in Section |3] In Section |4] we shift attention to the 
formation and evolution of debris fields. In this regard. Sec- 
tion U deals with mass loss from subhaloes and the contri- 
bution of so-called "backsplash" satellites. These subhaloes 
are found outside the virial region of the host at the present 
day, but their orbits took them inside the virial radius at ear- 
lier times. We show that they exist not o nly in cluster-sized 
haloes as found by e.g. I Gill et aL I (|2005h and iMoore et al.l 
(|2004l ). but also in our galactic halo. In Section[S]we examine 
how well satellite properties such as its original mass and the 
eccentricity of its orbit can be reproduced using the whole 
debris field of the corresponding satellite. Furthermore, we 
present in Section[6]a method to distinguish between trailing 
and leading arms in an automated fashion and we investi- 
gate properties of the streams such as their mass, shape and 
velocity dispersion. We pay particular interest to the energy 
distribution of debris, since it can be used to distinguish be- 
tween leading and trailing arm observationalljjf]. Finally, we 
summarise our results in Section [7] and comment on their 
significance in the context of future papers, which will try 
to establish the missing link between the properties of tidal 
streams and the host haloes only briefiy touched upon in 
this study. 



2 SIMULATIONS & ANALYSIS METHODS 
2.1 The Simulations 

Our analysis is based on a set of nine high-resolution cos- 
mological TV-body simulations that were used in the study 

For the purposes of our study, we treat substructure haloes and 
satellite galaxies as interchangeable, insofar as both suffer mass 
loss and the tidally stripped material is a dynamical tracer of the 
host potential. However, we note that the correspondence between 
dark matter substruct ures and lumino us satellite galaxies is not 
a straightforward one llGao et al.ll2004) . 

^ It is important to remember that our tidal streams consist of 
"particles" of dark matter, whereas an observer detects stellar 
streams. However, the physical processes involved in the forma- 
tion of both dark matter and stellar streams are identical, and so 
we can think of them interchangeably. 
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of lWarnick fc Knebd (|2006l ). These simulations focus on the 
formation and evolution of a sample of galaxy- and cluster- 
mass dark matter haloes forming in a spatially flat ACDM 
cosmology with f2o=0.3, f2A=0.7, f2i,=0.04, h — 0.7, and erg 
— 0.9. Each halo contains at least 10^ particles at z=0 and 
each simulation was run with an effective spatial resolution 

of S; O.OlTivir. 



Eight of the haloes - the cluster-mass systems, C1-C8 
- were simulated using the publicly availa ble adaptive mesh 
refinement code MLAPM (|Knebe et al.ll200fl ') . We first created 
a set of four independent initial conditions at redshift z = 45 
in a standard ACDM cosmology (Slo — 0.3, il\ = 0.7, Qt = 
0.04, = 0.7, o-g = 0.9). 512^ particles were placed in a 
box of side length 64/i~^ Mpc giving a mass resolution of 
nip — 1.6 X 10* Mq. For each of these initial conditions 
we iteratively collapsed eight adjacent particles to one par- 
ticle reducing our particle number to 128^ particles. These 
lower mass resolution initial conditions were then evolved 
until z = 0. At z = 0, eight clusters (labeled C1-C8) from 
our simulation suite were selected, with masses in the r ange 
1-3 X 10^* ft" ^ Mq. Then, as described bv lTormenI (|l997l ). for 
each cluster the particles within five times the virial radius 
were tracked back to their Lagrangian positions at the initial 
redshift (z = 45). Those particles were then regenerated to 
their original mass resolution and positions, with the next 
layer of surrounding large particles regenerated only to one 
level (i.e. 8 times the original mass resolution), and the re- 
maining particles were left 64 times more massive than the 
particles resident with the host cluster. This conservative 
criterion was selected in order to minimise contamination of 
the final high-resolution haloes with massive particles. 



The ninth (re-) simulation was performed using the 
same (technical) ap proach but with the ART code 
(jKravtsov et al.l Il997h . Moreover, this particular run (la- 
beled Gl) describes the formation of a Milky Way type dark 
matter halo in a box of side 20h ~^ Mpc. It is the sa me sim- 
ulation as "Box20" presented in lPrada et al.l ()2006[ ) and for 
more details we refer the reader to that study. The final 
object consists of about two million particles at a mass res- 
olution of 6 X 10^ h^^ Mq per particle and a spatial force 
resolution of 0.15/i"^ kpc has been reached. All nine simu- 
lations have the required temporal resolution (i.e. At « 200 
Myrs for the clusters and At ~ 20 Myrs for the galactic halo, 
respectively) to accurately follow the formation and evolu- 
tion of subhaloes within their respective hosts and hence 
are well suited for the study presented here. The high mass 
and force resolution of our simulations is sufficient to resolve 
an abundance of substructure (the "satellites" in our study) 
within the virial radius of each host. 



The physical properties of the haloes (based upon the 
halo finding technique explained below) are summarised in 
Table [T] alongside some additional indicators quantifying the 
mass assembly histories (cf. Section |3]). 



2.2 Halo Finding 

Both the haloes and their subhaloes are identified using 
AHtifl, an MP I parallelis e d mo dification of the MHljfl algorithm 
presented in lGill et al.l (|2004[ ). AHF utilises the adaptive grid 
hierarchy of MLAPM to locate haloes (subhaloes) as peaks in 
an adaptively smoothed density field. Local potential min- 
ima are computed for each peak and the set of particles that 
are gravitationally bound to the peak are returned. If the 
peak contains in excess of 20 particles, then it is considered 
a halo (subhalo) and retained for further analysis. 

For each halo (subhalo) we calculate a suite of canon- 
ical properties from particles within the virial/truncation 
radius. We define the virial radius i?vir as the point at which 
the density profile (measured in terms of the cosmological 
background density pb) drops below the virial overdensity 
Avir, i.e. M{< i?vir)/(47r7?vir/3) = Avirpb. FoUowiug con- 
vention, we assume the cosmology- and redshift-dependent 
definition of Avir; for a distinct (i.e. host) halo in a ACDM 
cosmology with the cosmological parameters that we have 
adopted, Avir = 340 at z = 0. This prescription is not 
appropriate for subhaloes in the dense environs of their 
host halo, where the local density exceeds AvirP6, and so 
the density profile will show a characteristic upturn at a 
radius R < -Rvir. In this case we use the radius at which the 
density profile shows this upturn to define the truncation 
radius for t he subhalo . Further details of this approach can 
be found in lGill et all (|2004l ). 



2.3 Subhalo Tracking 

In the present study we are primarily interested in the tem- 
poral evolution and disruption of subhaloes and hence a 
simple halo finder will not suffice. Therefore w e base our 
analys is on the halo tracking method outlined in et al.l 
(|2004l f1. At the formation time of the host halo we iden- 
tify all subhaloes in its immediate environment using AHF 
and track their subsequent evolution using the set of par- 
ticles identified at the initial time. Using the 25 innermost 
particles at time ti we calculate the new subhalo centre at 
the subsequent time t2. We then check which of the initial 
set of particles remain bound to the subhalo; all unbound 
particles constitute "debris". Particles are unbound if their 
speed exceeds the escape velocity by a factor of 1.5. This 
tolerance factor is introduced to allow for any inaccuracy 
arising from the assumption of spherical symmetry or un- 
certainties in the extent of the subhalo. All particles, bound 
as well as unbound ones, are tracked forward to the present 
day, which allows us to look at the field of debris particles 
for each subhalo. 

The tracker does not take into account accretion of 
mass, and so we may slightly underestimate subhalo masses. 
However, it is very unlikely that subhaloes will accrete ma- 
terial as they orbit within the host; rather, they will suffer 

^ AMIGA's-Halo-Finder; AHF can be downloaded from 
http : //www. aip . de/People/aknebe/AMIGA. AMIGA is the successor 
to MLAPM. 

MLAPM's-Halo-Findcr 
^ The tracking utility HaloTracker is part of the publicly avail- 
able AHF distribution. 
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continuous mass loss. In contrast to a simple halo finder, 
the tracker is able to detect "flickering" subhaloes that tem- 
porarily disappear in the dense inner regions of the host (e.g. 
Power & Knebe 2007). The tracker is in the exceptional po- 
sition of knowing that there should be a subhalo and thus 
can track it as it orbits within the densest innermost parts 
of the host fsee lGill et aLll2004l ). 

However, there is one draw-back with this otherwise 
useful feature: a subhalo may completely merge with the 
host, but its particles still appear to constitute a bound 
structure. Such 'fake' subhaloes appear to lie close to the 
centre of the corresponding host. Thus, we remove these ob- 
jects by neglecting every subhalo whose centre lies within 
the central 2.5% of the virial radius of the host and which 
remains within this radius for several consecutive time-steps. 



2.4 Orbit analysis 

For studying subhaloes and their properties, knowing their 
individual orbits around the host centre is important. We 
focus on subhalo orbits in the rest frame of the host (which 
encompasses all (bound) matter inside its virial region) . Be- 
cause of inherent uncertainties in the accuracy with which 
we can determine a halo's centre, the path of the host halo 
appears to follow a "zigzag" trajectory. Therefore, as an ini- 
tial step we smooth over the host's orbit before translating 
subhalo orbits into the host frame. This ensures that sub- 
halo orbits do not artificially "wiggle" when translated to 
the host frame because of uncertainties in the host's precise 
centre. 

In addition, we need to calculate orbital properties such 
as the number of orbits around the host centre and eccen- 
tricity. We base our eccentricity determination on finding 
the apo- and pericentre of a subhalo's orbit. For each apo- 
(va) and pericentre (r^) pair, the corresponding eccentricity 
can then be defined as: 



However, noisy density determinations (and the infiu- 
ence of other subhaloes) may introduce additional minima 
and maxima in the orbit, thus leading to artificially low ec- 
centricities mimicking nearly circular orbits. We get rid of 
these minima and maxima by applying a "sanity" test to 
each supposed apo-pericentre pair with an eccentricity be- 
low 0.1: we estimate the time the subhalo would need for 
one circular orbit around the host using its mean distance 
(r) and speed {v) based on the values at the corresponding 
apo- and pericentre: T = 2n{r}/ {v). If the time between the 
suspected apo- and pericentre is less than 30% of half of 
the estimated orbital time, then this minima/maxima pair 
is most probably not related to the motion around the host 
centre and thus is removed. In this way, we achieve a much 
more reliable determination of apo- and pericentres for each 
subhalo orbit. 

This also affects the number of orbits that we estimate, 
since it relies on the number of found apo- and pericentres. 
Additionally, partial orbits at the beginning and end of the 
studied time interval are added to the number of full orbits 
derived from the extrema of the orbit leading to non-integer 
orbits. 



3 THE HOST HALOES 

Some properties of the hosts used in t his paper have al- 
ready been presented in a previous study (|Warnick fc Knebg 
boOfih . However, for completeness we list relevant properties 
in Table [T] the meaning of listed quantities are explained in 
this section. Note that the shape measurern ents we present 
differ from those in IWarnick fc Knebd (|2006h because we em- 
ployed a new method for shape estimation (cf. Section IX}. 

3.1 Formation Times 

The standard definition of the formation time of a dark mat- 
ter h alo is based on the simple criterion of lLacev fc Cold (e.g. 
Il993l ). which determines formation time to be the time at 
which a halo's most massive progenitor first contains in ex- 
cess of half its present day mass. If we apply this "half-mass" 
definition to our sample, we find that our haloes formed be- 
tween 3 <^ Afform <, 8 Gyrs ago. 

However, we can define an alternative formation time 
that is based explicitly on a halo's merging history. If we 
determine the time of the halo's last major merger, then 
we can use the age of the Universe at this time to define 
the halo's formation time. When applying this criterion we 
note that the "half-mass" formation time actually provides 
a lower limit on formation time; for all host haloes no merger 
with a mass ratio greater than 1:3 occurred since the half- 
mass formation redshift. We decided to use the conservative 
age determination based upon the half-mass criterion. 

3.2 Mass Accretion 

As already introduced in IWarnick fc Kneb3 (|2006l ) we com- 
pute the dispersion of the fractional mass change rate 

1 "^^'^ f AKh /AM\\' 

where A'out is the number of available outputs from forma- 
tion Zform to redshift z = 0, AMi = \M{zi) - M{zi+i)\ the 
change in the mass of the host halo, and Ati is the respec- 
tive change in time. The mean growth rate for a given halo 
is calculated as follows 

AM \ 1 ""^^ AMi , . 

AtM I Nout - 1 ^ AUM, ■ ^ ' 

!=1 

A large dispersion (Jam/m is indicative of a violent formation 
history whereas a small dispersion refiects a more quiescent 
history. Our definition for formation time implies that the 
host halo's mass has doubled its mass between the formation 
time and z=0, and as such the inverse of the age of the host 
provides an estimate for the mean growth rate, which can 
be compared with the rate calculated from Eq. (|3|l. 

3.3 Shape of the Gravitational Potential 

A generic prediction of the CDM model is that dark mat- 
ter haloes are triaxial syste ms, that can be reasonably well 
approximated as ellipsoid s (jFrenk et al ] |l988l : IWarreii et al.l 
ll992l : lAllgood et al.ll200^ . etc.). 

Here we apply the method detailed in Appendix 1X1 to 
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Table 1. Properties of the host haloes in our simulations. Masses are measured in 10^*/i~^ Mq, velocities in km/s, distances in kpc, 
and the age in Gyr. Mvir is the virial mass, i?vir is the virial radius, Knax is the maximum of the rotation curve, -Rmax is the position 
of this maximum, age is the time since half of the present day mass was first assembled, 2form is the corresponding formation rodshift, 
Ci/5 = -Rvir/J^i/s is a measure of the concentration (where -R1/5 is the radius containing 20% of the virial mass), and A is the (classical) 
spin parameter. The following two columns quantify the mass assembly of the hosts (in Gyr~^ and Gyr~^, respectively) and we refer 
the reader to the corresponding Section l3.2l for details. The last three columns measure the the shape of the hosts based upon our novel 
technique introduced in Appendix |X] and applied to the potential of the hosts in Scction l3.3l as well as the dynamical state Q defined in 
Section |3^ 



our suite o f host haloes and d efine the triaxiahty in the usual 
way (e.g., iFranx et al1ll99lh 



(4) 



where the ellipsoidal axes a > b > c are closely related to 
the eigenvalues of the inertia tensor (of. Eq. ((BT)!). 

Not e that the values for T li sted in Table [1] differ from 
those in IWarnick fc Knebd (|2006l ). In the previous study we 
used all particles within the virial radius (cf. dashed circle 
in, for instance. Figure IC1|I whereas the present values are 
based on the novel method outlined in Appendix [X] which 
uses only particles within an ellipsoidal equipotential shell 
with a "mean" distance (cf. _Riso defined in Appendix IA2|) 
from the centre of fa i?vir (cf. solid lines in the upper panels 
of Figure [CTIl . 



3.4 Dynamical State 

We use the virial theorem in order to quantify the dynamical 
state of our host haloes. For a self-gravitating system of A'' 
(collisionless) particles it is given by 



1^1 
2 dt2 



■ 2i5kin + Ep 



(5) 



where I is the inertia tensor. The potential Spot and kinetic 
Skin energies are given by 



Skin = ^ y^^(miVi) 



pot 



IE 



(6) 
(7) 



where the summations are over all N particles. If we as- 
sume that the inertia tensor does not depend on time, then 
expression ([5]) reduces to the more familiar 



Q = 2Skin/|Spot| — 1 , 



(8) 



which we refer to as the "virial ratio". The summation to 
obtain both the kinetic and potential energy is performed 
over all particles inside an equipotential surfaces. We chose 



the one whose "mean" distance (cf. Riso defined in Ap- 
pendix 1X2)1 to the host's centre is i?vir, just as for the shape 
determination in Section [3.31 



3.5 Summary — Host Haloes 

All our host haloes appear to be relaxed and in dynamical 
equilibrium with a range of formation times and triaxiali- 
ties. Hence we have an interesting sample of hosts to study 
the formation and evolution of tidal debris from disrupting 
subhaloes orbiting within and about them. However, visual 
inspection of the temporal evolution of the hosts reveals that 
C8 is experiencing an ongoing triple merger, which explains 
why it has the largest deviation of Q from unity amongst all 
haloes in our sample. 



4 THE SUBHALOES 

Strong tidal forces induced by the host halo will lead to a 
substantial mass loss from the orbiting subhaloes. In this 
section we examine this mass loss in detail and investigate 
how it relates to the orbital parameters of the subhaloes. 
In particular we consider two distinct subhalo populations, 
"trapped" and "backsplash" subhaloes. Trapped subhaloes 
enter the virial radius of the host at an earlier time and 
their orbits are such that they remain within the virial ra- 
dius until the present day, being forever "trapped" in the 
potential of the host. Ba cksplash subhal oes were introduced 
earlier (cf. Introduction, iGill et al]|2005l ). Both trapped and 
backsplash subhaloes suffer mass loss, but it is instructive to 
ask how this lost mass contributes to the tidal debris field 
within the host halo. In particular, backsplash subhaloes 
lose a substantial fraction of their mass as they pa ss through 
the de nser environs of their host (on average 40%. lGiU et ai] 
(|2005l )). but how important is their contribution to the tidal 
debris field? These are interesting questions that will be ad- 
dressed in this section. 
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Figure 1. Minimum distance of a subhalo versus present dis- 
tance. Sublialoes in the lower right corner of the plot are 'back- 
splash galaxies': they have once been inside the host, but now are 
sitting outside the virial radius. Subhaloes in the upper right part 
have never been inside the host, but most of them are approach- 
ing the host halo (rf^— q = dmin; solid grey line), (compare with 
Figure 2 of Gill et al. 2005)) 



4.1 Trapped &: Backsplash Subhaloes 

As noted by both lGill etHI (j2005l ) and lMoore et alj (|2004h a 
prominent population of backsplash galaxies is predicted to 
exist in galaxy clusters. We now demonstrate that a back- 
splash population is also predicted to exist in galaxies. In 
Figure [T] we show for each subhalo its minimum distance 
with respect to the centre of the host halo as a function of 
its distance with respect to the centre at the present day; 
both are normalised with respect to the host's virial radius 
at redshift z — 0. 

Figure [1] reveals that backsplash subhaloes are also 
present in our galaxy mass halo. The backsplash popu- 
lation occupies the bottom right rectangle, with ciz=o ^ 



IR 



■vir,host,z — 



and dn 



^ li?vir 



For the clus- 



ters, about 50% of the subhaloes within a distance d G 
[fl!vir,host, 2 i?vir,host] ouce passcd through the virial sphere 
of their respective host. For the galactic halo Gl, this figure 
rises to 63%. We note that the concept of the backsplash 
subhalo can partly be ascribed to the difficulty in choosing 
a good truncation radius for the host halo: the virial radius 
is more or les s an arbitrary defini tion for the outer edge 
of a halo (cf. iDiemand et al.l [20061 ) and it has been shown 
by others that bound objects (in isolation) ma y as well ex- 
tend way beyond this formal point-of-reference (|Prada et al.l 
l2006l ). 

The numbers of backsplash subhaloes for each host are 
given in Table [2j along with the number of trapped sub- 
haloes. We note that the number of subhaloes is slightly 
larger than is given in IWarnick fc Knebd (|2006l ). refiecting 
differences in how subhaloes are identified. In the previous 
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Table 2. Number of subhaloes inside the virial region of the host 
(n), backsplash subhaloes between 1 and 2 virial radii (njjs) and 
the fraction of backsplash subhaloes in that region (with respect 
to all subhaloes between 1 and 2 virial radii). Only subhaloes 
containing at least 20 particles are counted. 
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Figure 2. Fractional mass loss per orbit since infall time as a 
function of infall eccentricity. Subhaloes were selected according 
to the criteria: ^ 1 orbit, npart ^ 20 particles today and they 
must reside inside the host today. 



study we used AHF, which identifies a "static" distribution of 
subhaloes in a single snapshot, w hereas in the pre sent study 
we use the halo tracker AHT (see iGill et al] |2004| . for more 
details). 

Interestingly, the fraction of backsplash subhaloes 
within one and two virial radii of the host varies between 
43% and 63% for all hosts, except for the youngest system 
C8. There we find a rather low fraction of only 19%. This 
host is not only the youngest, but also an ongoing triple 
merger providing the most 'unrelaxed' environment amongst 
our suite of hosts (cf. Q value in Table [1]). Thus maybe there 
are some potential backsplash subhaloes which simply did 
not yet have the time to leave their host again. 
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Figure 3. The cumulative mass function of tlie backsplasli galax- 
ies at redshift 2 = (heavy lines) and when using the mass at 
the infall time of each individual backsplash galaxy (faint lines). 

4.2 Mass Loss from Trapped Subhaloes 

Before investigating tidal debris fields it is necessary to 
quantify the actual mass loss. Are our subhalo systems actu- 
ally losing mass and how does this mass loss relate to orbital 
parameters such as infall eccentricity? To answer this ques- 
tion we show in Figure [2] that trapped subhaloes on more 
eccentric orbits experience in general larger mass loss. We 
plot the total mass loss since infall time normalised to the 
initial mass divided by the number of orbits as a function of 
orbital eccentricity. The orbital eccentricity is derived here 
from the first apo- and pericentre distance of the subhalo's 
orbit, after the subhalo entered the host region. This "infall 
eccentricity" might deviate strongly from the present value, 
since the eccentricity usually is evolving within time (see e . g . 
[Hashimoto et al.ll2003l : iGill et al.|[20oi ISales et al.ll2007l ) . 

While in general all of the trapped subhaloes lose about 
30% of their mass per orbit, some objects on radial orbits 
may lose up to 80% per orbit or even more. We like to note 
that only surviving haloes (^20 particles) enter these plots 
and completely disrupted ones were not considered, respec- 
tively. Furthermore, only objects with at least one orbit were 
taken into account for Figure [2] as we needed one orbit to 
get a reliable measure for the (infall) eccentricity and the 
number of orbits. 

4.3 Mass Loss from Backsplash Subhaloes 

The results of the previous subsection show that trapped 
subhaloes suffer mass loss whose rate is enhanced for sys- 
tems on more radial orbits. We now examine mass loss from 
backsplash subhaloes. 

In Figure [3] we show the cumulative mass function of 
all backsplash subhaloes at redshift z = Q (heavy lines) and 
compare with the mass function of these subhaloes at their 



Figure 4. Fractional mass loss from subhaloes since infall time 
versus the normalised minimum distance to the host. Only sur- 
viving subhaloes are plotted. Trapped subhaloes are presented by 
blue plus signs, backsplash galaxies are marked by orange dia- 
monds. 



time of infall onto their host (faint lines). The two sets of 
mass functions approximately preserve their shape and differ 
effectively only in horizontal offset - the z=Q mass function 
can be obtained by a shift to the left of the mass function 
at infall. Further inspection of the figure reveals that ev- 
ery backsplash galaxy loses on average half its mass as it 
travels within the virial radius of the host and out again, ir- 
respective of its original (infall) mass, i.e. the mass functions 
presented for each host halo in Figure |3] are self-similar. 

Although backsplash subhaloes may not count as "true" 
subhaloes of the host, they experience mass loss within the 
host and therefore contribute to the tidal debris fields that 
we study in Section [5] and Section |6] We verify this in Fig- 
ure |3J similar in spirit to Figure [51 where we plot the total 
mass loss from a subhalo over its full orbit since infall time 
as a function of minimum distance in terms of current host 
radius. We separate trapped and backsplash subhaloes by 
plotting the former with blue plus signs and the latter with 
orange diamonds. This figure demonstrates that subhaloes 
that have small pericentres with respect to the centre of the 
host experience greater mass loss, regardless of whether the 
subhalo is trapped or backsplash. We show also the average 
mass loss at a given minimum halocentric distance by the 
histograms in Figure [4] Interestingly, it is noticeable that 
backsplash subhaloes (orange) that grazed the outer regions 
of the host suffer greater mass loss than trapped subhaloes 
(blue) that reside in the same region. 

This figure also underlines our previous assertion that 
backsplash galaxies contribute to the overall tidal debris 
field. Table [3] quantifies this: at most 13% of the debris 
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Table 3. Number of particles inside the virial radius of the host (at z = 0): Up = total number, rip^sat = the number of particles in still 
existing subhaloes {^20 bound particles), rip debris = total number of particles in debris (including destroyed subhaloes and debris from 
backsplash galaxies), rip ^spl = debris particles from backsplash subhaloes alone and the last column provides the ratio of the latter two 
numbers. 



particles found within the host originate from backsplash 
subhaloefl. 

Table |3] gives the fraction of particles in subhaloes, de- 
bris and backsplash galaxies; our hosts' masses in general 
consist of about 5 to 10% (in mass) from surviving subhaloes 
while 10 or 20% of the mass is contained in debris particles. 
A mere 0.5% of the total host mass originates from debris 
of backsplash galaxies. We further make the following obser- 
vations about the subhalo and debris fractions of individual 
hosts. 

• We note that the subhalo fraction in C8 is high, re- 
flecting its involvement in a recent merger. Because a single 
massive subhalo persists (in fact, one of the most massive 
progenitor of the host) that has not yet merged, the subhalo 
mass fraction is close to 36%. 

• The debris fraction for host C6 amounts up to nearly 
25% because this host contains several massive subhaloes. 
We checked the initial subhalo mass function, i.e. the mass 
function of all trapped subhaloes at the formation time of 
the host, and indeed C6 contains initially more mass in sub- 
haloes than any of the other hosts. The subhalo destruction 
rate is similar for all hosts and so the high debris fraction 
reflects the number of massive subhaloes that were present 
within the host at its formation time, rather than reflecting 
enhanced disruption of subhaloes within the host. 

• In contrast, C7 is a host with a very low debris con- 
tent, although it contains an average number of subhaloes. 
A detailed investigation of its temporal evolution revealed 
that neither the number of subhaloes nor the rate of subhalo 



° We note a few instances where we find negative mass loss, i.e. 
subhaloes actually gained mass. This is due to marginal technical 
differences between the halo finder and the halo tracker: while the 
finder provides the tracker with the initial set of particles belong- 
ing to the (sub-)halo the tracker checks whether these particles 
are bound or unbound. While the finder uses local potential min- 
ima as the centre of the halo the tracker re-calculates the centre 
as the centre-of-mass of the innermost 25 particles. If there is now 
a marginal shift this can lead to differences in the evaluation of 
the potential and hence the odd particle can be considered un- 
bound by the tracker. However, because it is in fact bound it will 
"merge" with the halo again mimicking mass growth. We care- 
fully checked these instances and those objects appearing with 
negative mass growth are in fact subhaloes that did not suffer 
any mass loss at all. 



disruption varies significantly with time. This is consistent 
with the small measured debris fraction. 

It is instructive t o compare our f i nding s agai nst the re- 
sults o f, for example, iDiemand et al.l (|2006l ) and iGao et al] 
(|2004 ). The former investigated mass loss per orbit as a func- 
tion of pericentric distance (cf. their Fig. 17) and arrived at 
the same conclusion: the closer a subhalo gets to the host's 
centre, th e greater i t s mas s loss. A similar result was ob- 
tained bv lGao et al.l (|2004l ) who found that subhaloes with 
large pericentric distance retains a greater fraction of their 
initial (i.e. infall) mass (cf. Fig. 15). 

4.4 Summary — Subhaloes 

Our study of the subhalo populations in the hosts reveals 
that backsplash subhaloes are present in both galaxy and 
galaxy cluster haloes. Approximately half of the haloes be- 
tween 1 and 2 virial radii are actually backsplash subhaloes. 
Both trapped and backsplash subhaloes suffer mass loss as 
they orbit within the dense environs of the hosts - between 
30 to 80% of their mass per orbit. The rate of mass loss 
is enhanced for those subhaloes that follow highly eccen- 
tric orbits and for those subhaloes whose orbits take them 
to small halocentric distances. This is in good agreement 
with the findings of previous studies (e.g., iGao et al]|2004l : 
iDiemand et al.ll2006l ). 

This lost mass will contribute to the tidal debris field 
within the hosts. The fraction of a host's mass that comes 
from this lost mass ranges from a few percent (C7) to about 
25% (C6). This diversity reflects differences in the subhalo 
populations that were present at the formation time of the 
host, and differences in the accretion histories of the hosts. 
We note that the contribution of mass lost from backsplash 
subhaloes to the total debris inside the host halo region to- 
day amounts to at most 13%. However, we conclude that 
when observing the debris of subhaloes and searching for pos- 
sible corresponding subhalo galaxies, the region outside the 
virial radius of the host galaxy should also be considered. 



5 TIDAL STREAMS I: 

UNCOVERING SUBHALO PROPERTIES 

When a subhalo moves on its path around the centre of 
its host, tidal forces (mainly from the host) tear it apart. 
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usually generating the formation of a leading arm ahead 
of the subhalo and a trailing arm behind it. These arms 
are not always clearly distinguishable; sometimes they are 
broadened to such a degree that only a rather smooth sphere 
of unbound particles remains. How does the morphology of 
the debris field depend on the subhalo and the host? Can 
we learn anything about the host halo or the subhalo itself 
by studying only its debris field? 

At th is point it is instructiv e to consider the recent 
studies of iPenarrubia et al.l (|2006l ) and ISales et"aL I (|200'if ). 
IPefiarrubia et alT i 20061 ) investigated the formation and evo- 
lution of tidal streams in analytic time-evolving dark matter 
host halo potentials and showed that observable properties 
of streams at the present-day can only be used to recover 
the present-day mass distribution of the host - they cannot 
tell us anything about the host at earlier times. However, if 
the full six dimensional phase space information of a tidal 
stream is available, it is also possible to deter mine the to- 
tal ma ss loss suffered by the progenitor subhalo. ISales et al.l 
(|2007l ) showed that the survival of a subhalo in an external 
tidal field depends sensitively on its infall mass. One of our 
primary interests is the relation between the debris field and 
the total mass loss of the progenitor subhalo as well as its 
orbital eccentricity. 

In the following analysis we focus on subhaloes that are 
the progenitors for tidal debris fields that contain at least 25 
particles. We have taken care to ensure that low-mass sub- 
haloes are not artificially disrupted (because of finite time- 
stepping, force resolution, particle number) by verifying that 
the fraction of a subhalo's mass in the debris field does not 
correlate with the mass of the subhalo at infall. If numer- 
ics were playing a role, then we would expect lower mass 
systems to show preferentially higher mass fractions in the 
debris field. We further remade all plots showing analysis of 
the debris fields restricting the minimum number of particles 
in the progenitor satellite to A'^niin = 600. We could not find 
any substantial differences between results and hence ad- 
here to the original plots. We note that the plots obtained 
by using the more restrictive limit on the satellite's mass ap- 
pear like having randomly sampled the original data; that is, 
they show the same trends. However, we note that in some 
instances extreme outliers are removed. 



5.1 Orbital and Debris Plane 

Definition 

For an ideal spherical host potential the orbit of a subhalo 
(and therefore its debris field) is constrained to a plane (i.e. 
the orbital plane) because angular momentum is conserved. 
In non-spherical potentials, however, the orbit precesses out 
of the plane, and this lead s to a thickening of the debris 



field in this direction (e.g. Ilbata et al.l [2OO1I : tHclmi 20o3; 



Ijohnston et al ] |2005l : IPenarrubia et aLlbOOd ). In addition to 
this effect, we expect the debris fields of more massive sub- 
haloes and of subhaloes on highly eccentric orbits to show 
a greater spatial spread because they are more likely to lose 
more of their mass over a larger volume. 

In order to quantify the deviation of both the subhaloes' 
orbit and debris field from the orbital plane, we first deter- 
mine a "mean orbital plane", i.e. the best-fit plane to the 
orbit of the subhalo. Utilising all available orbital points for 
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Figure 5. Fitting a plane to the orbit of three sample subhaloes 
in host halo CI. The left panel shows the orbit projected to its 
best fit plane, the right panel shows a perpendicular view, directly 
looking at the edge of the plane (straight red line). An asterisk 
marks the position of the subhalo today (2 = 0). The host halo's 
centre is located at (0,0,0). 



each individual subhalo (since infall tim e), we applied the 
Levenberg-Marquardt algorithm (see e.g. IPress et al ] |200l 
Chapter 15.5) to minimise the distances of the points (with- 
out loss of generality) in ^-direction. When using the general 
equation of a plane E : ax -\- by + cz + d = (with the nor- 
mal vector n = (a, &, c)), we obtain for a given x and y the 
expected z- value: 



z 

with 



ax -\- I3y -\- 5 
a 



5 = - 
c/0 



(9) 
(10) 

(11) 
(12) 



Employing the Levenberg-Marquardt algorithm then deliv- 
ers the coefficients a, jS and 5 describing our best fit plane. 

Examples of our fitting results are given in Figure (5) 
They illustrate how we can fit a plane to even a precessing 
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orbit (cf. bottom panel). We have inspected visually many of 
our fits to orbital planes and we find that the "mean orbital 
planes" provide an adequate approximation 0. 



Scatter about the orbit plane 

Given the best-fit orbital plane, we now have a means to 
quantify the deviation of the orbit itself from this 'mean' 
plane as well as the scatter of the debris particles about 
this plane. Using the orbital plane given by the equation 
E : ax + by + cz + d = 0, we determine the perpendicular 
distance Ai of each orbit point/debris particle to the plane 
via the (general) equation 



5,. = 



axi + byt + czi + d 



+ fo2 + c2 



(13) 



The mean deviation of the orbit Aorbit is determined using 



^orbit 



\ " c2 



3 ^ 

i=0 



St 



(14) 



The factor (n — 3) corresponds to the number of degrees of 
freedom, i.e. 'number of data points (n)' — 'number of free 
parameters (3)'. 

For the scatter of the debris, we use an analogous for- 
mula for the mean deviation: 



1 



^debris 



(15) 



Figure |6] shows the resulting debris scatter normalised 
to the deviation of the orbit itself for all chosen subhaloes 
in each simulation as a function of the infall eccentricity 
and the infall mass. Only subhaloes that have undergone at 
least one complete orbit inside the host halo are shown, since 
only then a reliable (infall) eccentricity determination based 
upon (first) apo- and pericentre of the orbit is possible. Dis- 
rupted subhaloes, or rather their debris, are also included 
in these plots. Subhaloes that once resided inside the host 
but are today found outside the host's virial radius, i.e. the 
backsplash galaxies (see Section l4.ip . are marked by grey 
asterisks. Interestingly, these are often the points deviating 
most from our expectation: we expect to observe a corre- 
lation of deviation with (infall) eccentricity, i.e. the larger 
the orbital eccentricity the greater the deviation of the tidal 
debris from the orbit. 

Only a marginal trend is apparent that is most pro- 
nounced for the oldest systems (i.e. CI and Gl): CI, on 
the one hand, is one of those with a prolate shape, deviat- 
ing more strongly from sphericity than most others. Thus, 
the large deviation of debris from the orbit plane could 
be a result of the flattened host halo as suggested by con- 
trolled experiments of disrup ting subhaloes in flattened an- 



alytical host potentia l s (e.g. Ibata et al 



Ijohnston et al.l l2005l : iPenarrubia et al 



.20031 : lHelmill2004 
20061 ). Gl, on the 



other hand, contains the most subhaloes, thus encounters of 
debris streams with other sub-clumps are very likely, caus- 
ing a heating of the stream which spreads the debris even 



^ The interested reader can look at Figure [TTI in Section |61 which 
shows examples of not only the subhalo orbit but also the debris 
field in projections of the best-fit orbital planes. 



more (iMoore et al.lll999l: llbata et al.|[20o3 : IJohnston et all 
l2002l : IPeiiarrubia et aLlbOOd ). 

The right panel of Figure [6] suggests that the trend is 
significantly more pronounced when plotting the deviation 
from the orbit plane against (infall) mass of the subhalo. 
Massive subhaloes always lead to large deviations of the de- 
bris about the orbital plane, whereas the debris of low mass 
subhaloes shows all kinds of deviations. We cannot construct 
a unique mapping of "debris deviation" onto "progenitor 
mass" by this method alone, but we no te that it agrees with 
the findings of lPefiarrubia et~aL l (|2006t ). who claim that it is 
possible to infer the amount of mass loss from the stream's 
progenitor using the full phase-space distribution of a tidal 
stream. We will return to this point later when studying the 
energies of stream particles. 



Scatter about the debris plane 

An observer will not be able to recover the orbital path 
of a subhalo galaxy (without, for instance, semi-analytical 
modeling of its trajectory) and thus no information about 
the orbital plane will be directly available to him (or her). 
However, he (or she) could either define a plane based on 
the present subhalo's position and velocity (i.e. the plane 
defined by the angular momentum vector) or fit a plane to 
the detected debris field and then determine the deviation 
of the debris field from this best-fit plane. We have adopted 
the latter approach and the result can be viewed in Figure [7] 
where we show the scatter of the debris particles about the 
best-fit debris plane versus infall eccentricity and mass. We 
note that the deviations of debris from a plane defined in 
this way are smaller than the deviations from the orbital 
plane. Again, only a marginal trend of increasing deviation 
for eccentric orbits is apparent. However, we now observe a 
stronger correlation between infall mass and debris devia- 
tion. Thus, by measuring the deviation of debris from the 
debris plane itself, it is possible for an observer to predict the 
most probable mass of the subhalo; in contrast, recovering 
the subhalo's eccentricity will be more ambiguous. 

The question arises as to why the correlation between 
(infall) mass and the debris variation is stronger in the case 
of the debris plane when compared to the orbital plane. Fur- 
thermore, we do not observe an "improvement" in the cor- 
relation with orbital (infall) eccentricity. We attribute this 
to the fact that the stream particles follow naturally the 
debris plane and hence show greater deviation from the or- 
bital plane (that does not necessarily coincide with the de- 
bris plane); the debris plane is defined such that it minimises 
the scatter of the stream particles about that plane. There- 
fore, the correlation already apparent in Figure |6] strength- 
ens when reducing the scatter in the deviation as is done 
when considering the debris plane rather than the orbital 
plane. 

One may further argue that the destruction of subhaloes 
depends on their (infall) concentration and t hat more mas- 
sive s ystems have lower concentrations (e.g., iBuUock et al.l 
l200lD . We confirm the expected scaling between mass and 
concentration, although we find no correlation between (to- 
tal) mass loss and (infall) concentration. Therefore, the 
trends seen in Figure [6] and [7] do not represent simple man- 
ifestations of the mass-concentration relation. 
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Figure 6. The deviation of debris particles from the best fitting orbital plane (with respect to the deviation of the orbit itself) as a 
function of infall eccentricity (left panel) and infall mass (right panel) for subhaloes with at least one orbit. Backsplash subhaloes are 
marked with grey asterisks. The median values of the distribution (trapped and backsplash galaxies together) are given in the lower right 
corner of each panel. 
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Figure 7. The deviation of debris particles from the best fitting debris plane (normalised to the virial radius of the respective host) as 
a function of infall eccentricity (left panel) and infall mass (right panel) for subhaloes with at least one orbit. Backsplash subhaloes are 
marked again with grey asterisks and the median values of the distributions are shown in the lower right corners. 
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Figure 8. The tube radius as a function of infall eccentricity (left panel) and infall mass (right panel), normalised to the virial radius 
of the host today. Only subhaloes with at least one orbit and with a well-defined tube radius (# particles inside tube = 68.3% ± 2%) 
are plotted. Grey asterisks correspond to backsplash subhaloes again. Numbers in the lower right corner of each panel show the median 
values for both, trapped subhaloes and backsplash subhaloes together. 



5.2 Tube Radius 

Even though the deviation from the orbital and debris plane 
may provide a means to infer the (infall) mass of a dis- 
rupting subhalo galaxy, we do not know whether the de- 
bris field is confined to the or bit or becomes c omple tely 
detached from it. The work by iMontuori et al.l (|2006h on 
tidal streams of globular clusters suggests that, in general, 
streams are good tracers of th e orbit of their pr ogenitor, at 
least on large scales. However. IChoi et al.l |20o3) argue that 
the mass of the progenitor has a great influence: a massive 
subhalo gravitationally attracts the debris field, pulling it 
back and therefore causing a stronger bending of the arms. 
Thus, the streams should deviate strongly from the orbital 
path, the larger the mass of the progenitor subhalo is. 

Furthermore, a large difference between the orbit and 
tidal streams may also indicate that the debris is just spread- 
ing over a larger area. This could then again be related to the 
infall eccentricity: subhaloes on spherical orbits will proba- 
bly keep their streams quite close to the orbit, whereas sub- 
haloes on radial orbits could have their debris spread to a 
much wider degree. 

To quantify how debris relates to the orbit we intro- 
duce a "tube analysis" - we generate a tube about the or- 
bital path of the subhalo, adjusting its radius rtubc such that 
1(7 (68.3%) of all debris particles are enclosed by the tube. 
Larger tube radii would then indicate that debris particles 
are deviating more strongly from the orbit of the subhalo. A 
detailed explanation of the "tube analysis" method and how 
to overcome its complications is presented in Appendix ID] 

The resulting tube radii (normalised to the host's virial 
radius) as a function of eccentricity and mass are plotted 
in Figure [S] Subhaloes are ignored, if the tube radius did 



not converge satisfactorily (more than 2% deviation from 
the desired la fraction of particles, mostly due to a small 
number of unbound particles). Again, only subhaloes with 
at least one full orbit are taken into account. 

The expected trend of increasing tube radii with mass 
and eccentricity is a weak one at best and we conclude that 
such a tube analysis is less helpful in revealing subhalo prop- 
erties than investigating the spread perpendicular to the de- 
bris plane. However, it nevertheless shows that the devia- 
tion of debris from the orbital path can be quite large. We 
therefore conclude that debris is not a well defined tracer of 
the orbital path of the progenitor subhalo as suggested by 
IChoi et all (|2007l ). 

We note also that in addition to studying how tube radii 
relate to subhalo properties, we examined how median tube 
radii related to the shape of the host. However, we were 
unable to identify any correlation. 



5.3 Radial Velocity 

So far we have utilised the 3D position (and velocity) in- 
formation available to us, but is interesting to ask how a 
subhalo stream would appear to an observer. An observer 
can, in principle, measure radial velocities of material in the 
tidal stream of a subhalo as well as its position on the sky. 
Let us imagine an observer positioned at the centre of a 
host halo, measuring the radial velocity of the (debris) par- 
ticles of a certain subhalo in our simulations. We fix our 
coordinate system such that the plane best fitting the de- 
bris lies at zero latitude (i.e. zero height). A longitude of 
shall correspond to the (remnant) subhalo's position today. 
This enables us to plot an 'observed' radial velocity distri- 
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Figure 9. The radial velocity of subiialo #32 in host CI (nor- 
malised to the maximum of the host's rotation curve) versus lon- 
gitude for an observer sitting in the centre of the host halo. Grey 
dots are bound particles, the black line shows the orbital path 
of the subhalo. Leading arm particles are coloured blue, trailing 
arm particles are shown in red with the brightness corresponding 
to the time since the particle became unbound (see Section |6] for 
more details). Note that we periodically wrapped the longitude 
extending the presented range from -180 to 360. 

bution versus longitude for a subhalo, as shown in Figure [Q] 
Note that only debris particles lost after the infall of the 
subhalo into the host are considered; see Section |6] for more 
details. From this figure it is apparent that the arms are 
following a sine-like curve, just like the orbit (black line). 
The trailing arrrF°l (red) sits neatly on top of the orbit path 
in this example and the leading arm (blue) seems like a 
perfect continuation of the orbit to future times. This par- 
ticular subhalo corresponds to one of the cases in the tube 
analysis that had a small tube radius. Its particulars are 
N^Mt = 2.2, eccinfaii = 0.458, Afi„faii = 8, M,=o = 1.6, 

Mtrailing arm = 3.7, Mlcading arm = 2.7 (whcrC maSSCS are 

measured in 10"/i"^ M©). 

Such plots contain a remarkable wealth of information 
about the subhalo and its orbit. The maximum/minimum 
radial velocity are a direct measure of the eccentricity of 
the orbit: for circular orbits we would expect a constant line 
at Vi = 0; the more eccentric an orbit is, the larger are 
its minimum/maximum radial velocity (in absolute values). 
Furthermore they hint at the number of orbits a subhalo has 
already undertaken: each period which the sine-like curve of 
the debris describes corresponds to one half (debris) orbit. 
The phase-shift of successive sine-curves in these plots could 
also be used to infer the precession of tidal debris in the de- 
bris plane, which may be correlated with the triaxiality of 
the host halo. Also, the spread of each arm in radial velocity 
could also tell us something about the mass of the subhalo's 
progenitor: more massive subhaloes will have a broader dis- 
trib ution of debris, in spatial terms as well as radial velocity 
(e.g. lChoi et"aLll2007l ). 

However, for the moment we concentrate on the total 
spread of the radial velocity 

A max mill / 1 r'\ 

Vr = — Vj. , (ib) 

and its relation to infall eccentricity and progenitor mass 
again. Figure \W\ shows the difference between smallest and 

^'^ Details on how to separate trailing and leading arm will follow 
in Section |6] 



largest radial velocity measured for the corresponding sub- 
halo's debris as a function of infall eccentricity (left panel) 
and initial subhalo mass (right panel). Indeed we observe the 
expected trend, although the scatter is very large. We can 
also see that subhaloes with more than one orbit (darker plus 
signs in Figures [TO} show a stronger signal: the more orbits 
undergone, the greater the probability that the difference in 
maximum and minimum radial velocity is determined by the 
orbital eccentricity rather than random scattering of the de- 
bris. But we also need to bear in mind that the debris field 
originating from subhaloes that have undergone a greater 
number of orbits will tend to be more massive, and there- 
fore we would expect to make a more precise determination 
of Avr- 

Using a more elaborate method for estimating the 
spread in radial velocity, such as using the mean values 
from the ten largest and smallest velocities or binning the 
distribution and taking the largest and smallest bin with a 
certain minimum number of particles, does not improve our 
results. 

We now turn to the correlation between the spread in 
radial velocity and the progenitor's infall mass (right panel 
of Figure llOp : we do not expect to observe a strong cor- 
relation, because only the width of the streams themselves 
should depend on the mass, whereas the total spread corre- 
lates more strongly with orbital parameters. However, Fig- 
ure [To] reveals quite a pronounced mass dependence: more 
massive haloes show a larger variety in radial velocities. Yet 
this correlation vanishes when looking at the radial velocity 
spread of the orbit itself (not shown here though). Thus, 
the signal found in Figure [TOl suggests that the debris does 
not merely follow the orbit but spreads out to some angu- 
lar extent that depends on its mass (compare also with the 
example radial velocity plot in Figure [Q]). Again this shows 
that the stream itself is more sensitive to the mass of its 
progenitor subhalo than the orbital properties, as already 
suggested by the results of the previous subsections. 

5.4 Summary — Uncovering Subhalo Properties 

We set out to recover properties of the progenitor subhalo 
from its tidal debris field alone. Therefore we focused upon 
correlations between the tidal debris field and the progenitor 
subhalo's infall eccentricity and mass. 

Although the expected correlation with infall eccentric- 
ity is weak, the statistics that we have presented may nev- 
ertheless help to recover the orbit of an observed subhalo. 
Both methods (scatter about the debris plane and tube anal- 
ysis) rely on measuring the debris field, thus allowing for a 
determination of the subhalo's eccentricity. We therefore en- 
vision that a combination of both the debris scatter about 
the best fit debris plane and the tube analysis applied to 
existing stream data may help to refine models for the re- 
construction of orbital parameters. 

Relating the tidal streams with the original mass of the 
subhalo, we observed a noticeable correlation when using 
the deviation of the debris from the mean debris plane. This 
opens the opportunity to infer the original mass of a (cur- 
rently) disrupting su bhalo such as , for in stance, the Sagittar- 
ius dwarf spheroidal llbata et al.l (|l994 ) . Using the relation 
between mass loss and eccentricity as presented in Figure [2j 
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Figure 10. The spread in radial velocity of unbound particles (normalised to the maximum velocity of the rotation curve of the 
corresponding host) as a function of infall eccentricity (left panel) and normalised initial subhalo mass (right panel). The grey crosses 
correspond to subhaloes with A^orbit < 3 while the black crosses represent subhaloes with A^orbit 3. 



it may be possible to reconstruct the eccentricity of the orbit 
indirectly as soon as we know the original and present mass 
of the corresponding subhalo. 

In addition, we presented an analysis of radial velocity 
plots as they might appear for an observer in the centre of 
a host. The spread in the radial velocity of the orbit itself 
is obviously strongly correlated with the (infall) eccentricity 
(not shown here), but it is less pronounced for the spread 
in debris (cf . Figure llOfl . It nevertheless shows a correlation 
with the mass of the progenitor, which is attributed to the 
fact that the debris of massive subhaloes actually spreads to 
a greater degree than for low-mass objects. 



6 TIDAL STREAMS II: 

SEPARATING LEADING & TRAILING ARMS 

So far the debris particles have been characterised by "bulk" 
properties, such as deviation from the orbital plane. No dis- 
tinction has been made between leading and trailing arm or 
the time since the debris particles became unbound. In this 
section we investigate differences in the properties of both 
the leading and trailing arms, which makes it necessary to 
tackle the rather challenging task of assigning particles to 
the trailing or leading debris fields. 

6.1 Age 

We begin by defining the age of stream particles, in order 
to understand the build-up of the stream. We classify a par- 
ticle as belonging to the stream at the moment at which it 
becomes unbound from its progenitor subhalo, which occurs 
when their velocity exceeds 1.5 times the escape velocity of 



the subhalo (see section . This age information is used 
to colour-code some of the following plots. 



6.2 Membership to leading/trailing arm 

We require not only the age of stream particles but also 
a means to determine which of the arms, i.e. leading or 
trailing, the particle belongs to. For subhaloes with clearly 
distinguishable arms, one can separate leading and trailing 
stream easily by eye; however, an automated method that 
can objectively classify particles is both desirable and nec- 
essary, especially when analysing literally hundreds of sub- 
haloes. We have developed such a method, which we now 
present. 

For each snapshot, we mark particles leaving the sub- 
halo in the "forward" direction as leading and those leaving 
in "backward" direction as trailing arm particles, where we 
define "forward" and "backward" as follows. If a particle 
leaves the subhalo, then - in a simplified picture - the sub- 
halo has lost its infiuence on this particle; the host now plays 
the most important role and hence all subhaloes shall be 
neglected. In addition, the velocity of a recently unbound 
particle will still be close to the velocity of its progenitor 
subhalo. In the absence of any forces the particle would fol- 
low trajectory defined by the direction of the velocity vector, 
but it feels the gravitational pull of the host and therefore 
its direction is defined relative to the vector 

^^^^_GA4 ,(r)r 
r V r 

If a particle leaves the subhalo in this direction, it is tagged 
leading, otherwise trailing. For a subhalo moving at high 
velocity, a stream will form mostly along the orbit (i;^-term 
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Figure 11. Classifying stream particles according to age and leading/trailing arm. Shown are the particles of a representative sample 
of subhalo galaxies orbiting in host CI projected to their orbital planes (left panel) and as an edge-on view (right panel). The black 
line represents the path of the orbit. Grey particles are still bound to the subhalo, red particles are marking the trailing, blue ones the 
leading arm. The colours vary from dark to bright, indicating the ago of the respective stream. 



prevails), whereas the stream of a slow subhalo is more at- 
tracted by the host centre (see e.g. iMontuori et al.ll2006l ). 

We note that using either the direction of the subhalo's 
velocity (v) or the direction to the host (r) alone when dis- 
tinguishing between leading and trailing arm particles does 
not give satisfactory results. Streams cannot be separated 
in this way because particles are not leaving simply in the 
direction of r or v. 

Our leading/trailing classification method works well, 
as indicated by the example subhaloes presented in Fig- 
ure [11] (see also the colors in the radial velocity plot. Fig- 
ure [9] where the same classification procedure was used). It 
fails, however, for those subhaloes which have just recently 
fallen into the host. Here the subhaloes may have already 
lost mass before entering the host, due to interactions with 
adjacent subhaloes or other overdensities. Particles lost in 
such a way cannot be described by our (nevertheless simple) 
picture since the influence of other subhaloes is too large 
compared with the host halo's influence. We therefore re- 
strict the remaining analysis of debris to particles lost after 
their (progenitor) subhalo entered the host. 

We note that Figure [TT] provides further evidence that 



our procedure for fitting orbits to a plane as outlined in 
Section rs . 1 1 works well for a variety of orbits. 

6.3 Masses in streams 

We are now in the position to explore differences in trailing 
and leading arm in detail and start with explicitly calculat- 
ing the ratio of the mass in both arms. Figure [T2l shows the 
ratio of the number of particles in the leading and trailing 
arm versus the infall eccentricity and mass of each subhalo. 
Again only subhaloes inside the host today and with at least 
one complete orbit (inside the host) are taken into account. 
Further, only subhaloes with at least 25 particles per arm 
are considered 

There appear to be slightly more particles in the 
leading than in the trailing arm in each of our cluster 
sized haloes while the reverse is true for the galactic halo 
Gl. This is interesting because of the recent claim of 

Interestingly, when including subhaloes with less massive 
arms, the number of points below 1 increases (more subhaloes 
with massive loading arm)! 
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Figure 12. Ratio of tlie number of leading and trailing arm particles for subhalocs inside the host with at least one orbit, plotted versus 
the infall eccentricity (left panel) and infall mass (right panel) . The number in the lower left corner of each panel gives the ratio of the 
number of subhaloes with a value above and below 1 . 



iKesden fc KamionkowskH (|2006l ). who argue that if dark 
matter is coupled to dark energy, there must exist another 
force acting on dark matter particles besides gravitation. 
This additional force could manifest itself in non-equilibrium 
systems such as the tidal tails of subhaloes, leading to an 
imbalance in the mass of the leading and trailing arms. 
According to a simplified picture of tail formation in cos- 
mologies without such coupling, both arms should contain 
the same amount of mass. Even though our simulations do 
not couple dark matter to dark energy, we observe a trend 
for an imbalance in the masses of leading and trailing arm. 
We therefore conclude that the proposition put forward by 
iKesden fc KamionkowskH (|2006l ) cannot hold in general. 

In addition, we do not observe any tendency for the 
ratio of masses in leading and trailing arms to correlate with 
infall eccentricity, but we do note a trend for more massive 
sub haloes to l o se mo re of their mass to the trailing arm. 
As IChoi et"al] (|2007l ) have pointed out, the gravitational 
attraction of the subhalo on the tidal arms may cause an 
asymmetry in leading and trailing arm. The strength of this 
effect will increase with the mass of the subhalo, which can 
explain the mass dependence we observe. However, it is not 
clear why this ratio reverses for the galactic halo, although 
we note that the reversal is evident only for the least massive 
subhaloes. 



6.4 Velocity dispersion in tidal streams 

The velocity dispersion of streams is another property that 
can reveal important insights to the history of the subhalo 
and ch aracteristics of the host. For example. I Johnston et al.l 
1 19961 ) note that a high velocity dispersion in the tail is 
suggestive of a massive progenitor as well as a high internal 
velocity dispersion in the original subhalo. However, the dis- 
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Figure 13. The velocity dispersion of trailing (black) and leading 
(dark grey) arm, normalised to the velocity dispersion of the host, 
as a function of the number of orbits. Since 'partial' orbits at the 
beginning and end of the considered time interval are also taken 
into account, the number of orbits is not limited to integers. 



persion can increase further if the stream orbits in a "lumpy" 
host halo: encounters between the stream and sub-clumps 
will lead to dynamical heating of the stream. Furthermore, if 
the shape of the host deviates strongly from spherical sym- 
metry, the subhalo and thus also its debris particles will 
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Figure 15. Ratio of the velocity dispersion in leading and trailing arm for subhaloes in the host with at least one orbit as a function 
of infall eccentricity (left panel) and infall mass (right panel). 
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Figure 14. The velocity dispersion of trailing (black) and leading 
(dark grey) arm as a function of the deviation of debris from the 
best fit debris plane (compare with Figure [7}. 



This figure is accompanied by Figure [14] where we show the 
relation between the velocity dispersion of debris and its 
scatter about the best fit debris plane (cf. Section I5.1|l . It 
underlines the trend that a larger deviation of the debris 
from the best fit plane correlates with a higher velocity 
dispersion and vice versa. Viewed in conju nction with Fig- 
urej6] (right panel) we recover the findings of lJohnston et ahl 
(jl996l l that a larger progenitor automatically leads to a 
higher velocity dispersion. 

Are there differences in the velocity dispersions of the 
leading and trailing arms, induced by the host halo? In Fig- 
ure [15] we show the ratio of the velocity dispersions in both 
tidal arms as a function of the subhalo's infall eccentricity 
and mass again. We observe that the mean ratio is always 
of order unity with a roughly constant scatter when plot- 
ting against infall eccentricity, but the scatter decreases for 
initially more massive subhaloes (right panel of Figure I15|) . 
There appears to be only marginal differences between the 
leading and trailing arm with respects to the internal ve- 
locity dispersion. As they are spatially separated entities we 
therefore suspect interactions with sub-clumps in the host 
halo to have little infiuence unless both arms accidently en- 
counter the same number of subhalo. It is more likely for 
the host to play the dominant role in the determination of 
the velocity dispersion of these two detached debris field. 



start to precess, leading to yet further spread of the debris, 
increasing its spatial dispersion. 

We start by exploring the relation between the number 
of orbits and the velocity dispersion in both the leading and 
trailing debris arm. The result can be viewed in Figure 1131 
where we find a general increase in the velocity dispersion 
with the number of orbits - the greater the number of 
orbits, the higher the velocity dispersion of the debris field. 



6.5 Energy of tidal debris 

We have learned that the masses of leading and trailing arm 
can differ, but their velocity dispersions are similar. While 
this may make it difficult to separate leading and trailing 
arms using raw velocity dispersions, we wish to understand 
whether we can use kinematic information to separate parti- 
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cles (or stars) into leading and trailing arms, as an observer 
might wish to do. 

To this end, we present a method that allows classifi- 
cation into leading and trailing arms based upon the en- 
ergy an d angular momen tum of tidal debris. Inspired by the 
work of IJohnstonI (| 19981 ) who found a bi-modal energy dis- 
tribution representative of leading (decrease in energy) and 
trailing arm (increase in energy) we calculate fo r each de- 
bris p article its change in t otal energy (also see McGlvniil 
Il99d : iJohnston erall 1 19991 . l200ll : iPefiarrubia" et al. I l2006t ). 
When combining the resulting distribution with our knowl- 
edge about the association of p articles with lead ing/trailing 
arms we recover the findings of Ijohnstoni (|l998l ) that there 
is a clear separation between the two arms in all our 'live' 
host haloes - at least for the majority of subhaloes. In what 
follows, we focus on the properties of streams at z—0. 



Energy of debris particles 

The energy of a particle is simply the sum of its potential 
and kinetic energy. While it is straightforward to calculate 
the kinetic energy, the potential energy is always somewhat 
more complicated. We are left with the option to either use 
the AHF analysis of the simulation which leaves us with (more 
or less) spherically cut dark matter haloes or go back to 
the original simulation and derive the full potential used to 
integrate the equations of motion. As the former relies on 
assumptions that are not necessarily true for all our sys- 
tems (i.e. we know quite well that our cosmological dark 
matter haloes are not spherical and better modeled by a tri- 
axial density distribution) we chose to follow to the latter 
approach. This provides us with the most accurate measure 
for the potential energy we can get. We also mention that 
we estimate both the potential and kinetic energy in the rest 
frame of the underlying host halo. We understand that this 
approach works against any hypothetical observer, but we 
wish to maximise the possibility that there is any chance at 
all to separate particles into leading and trailing arm. 



Energy scale 

According to IJohnstonl (| 19981 ) the changes in the total en- 
ergy are of the order of the orbital energy of the subhalo 
itself. In order to scale the variation in the particles' ener- 
gies in a similar fashion we define this scale to be 



AE E 



e — Epot + fkin 



(18) 



where epot and ekin are the specific potential and kinetic en- 
ergies of a test particle at the position and with the velocity 
of the (remnant) subhalo at redshift z — 0. Note that for tpot 
we interpolate the potential as returned by the simulation 
code AMIGA to the test particles' position. 



Energy distribution 

In order to be able to make predictions for an observer with 
access to the full 6D phase-space information for debris par- 
ticles/stars and a realistic model for the host halo (in order 
to derive potential energies), we restrict the following anal- 
ysis to quantities measured at z = 0. We define the relative 
energy change of a debris particle to be 



(19) 

c c 

where E represents the particles total energy and e the en- 
ergy scale as defined above. 

The resulting distribution for the same sample of sub- 
haloes as already presented in Figure [TT] can be viewed in 
Figure flBl This figure actually contains a substantial amount 
of information beside the energy distribution and hence re- 
quires careful explanation. 

First of all. Figure [T6l show s the distr i bution of q as de- 
fined by Eq. (fT^ . According to Ijohnstoni (| 19981 ') we should 
find a bi-modal distribution, at least for those subhaloes that 
are not yet completely disrupted. We therefore added the 
fraction of mass in the stream of the total mass to the legend 
(i.e. nu/ntot, where Uu is the number of unbound particles 
and Titot the total number of particles left in the subhalo and 
Uu together). To further highlight non-disrupted subhaloes 
the legend box in the upper right corner is shaded in grey for 
objects where less than 80% of all particles are forming the 
debris field (iiu/ntot < 0.8). As we are in the unique position 
to separate leading and trailing arm debris using the method 
described in Section 16.21 we not only present the combined 
energy distribution (black histograms) but also the individ- 
ual distributions for trailing (red histograms) and leading 
arm (blue histograms). The colour-coded dots above the his- 
tograms are the relative energies q of each individual debris 
particle where the saturation of the colour corresponds to 
stream age as defined in Section [6.21 too. 

While many of the energy distributions are clearly bi- 
modal (even for disrupted subhaloes), we notice that the 
dip in between the peaks is not necessarily located at g = 0, 
which just means that our choice of the energy scale is not 
good enough. Nevertheless, we claim that (i) this distribu- 
tion can be obtained observationally once the full 6D phase- 
space information (for at least a substantial amount of debris 
particles/stars) along with a sufficient model for the underly- 
ing host potential is available and (ii) subsequently be used 
to separate leading from trailing arm even if the debris is 
wrapped around the host and the progenitor subhalo had 
a number of orbits, respectively. One may claim that our 
model for the host potential is out of reach for an observer 
as we are using the exact potential values as returned by the 
simulation code, but we also applied the less sophisticated 
approach of using a simple spherically symmetric potential 
and were unable to detect significant differences in the re- 
spective energy distributions. 

Another noticeable feature in Figure [TB] are the 
"swings" in the particle distribution: these most probably 
correspond to the orbital motion of the subhalo. If we were 
scaling the particle's energies with the corresponding en- 
ergy scale when the particle became unbound or with the 
corres ponding pericentre energy (as suggested bv Ijohnstoni 
(|1998D '). these "swings" might vanish and we could possibly 
achieve a "cleaner" distribution. However, this energy scale 
would not be directly accessible to an observer. 



Separating leading and trailing arm via energy 

In order to check the credibility of our proposed method to 
classify a particle/star by its (relative) energy and assign 
it to either the leading or trailing arm we need to check 
the statistical significance of this technique. To this end we 
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Figure 16. The distribution of energy change of debris particles for a random sample of subhaloes in the host CI. Please refer to the 
main text for details. 
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Figure 17. Reproducibility of leading and trailing arm energy 
distribution by decomposing a double-Gauss-fit: The goodness 
for trailing (black) and loading (grey) arm is plotted versus the 
fraction of debris particles. 



CI 

-1- 
-1- 


C2 

-1- 

-Kf 

+ + + + ^ 


C3 

+ 

+ + + 

+ + ^ * 


C4 

-1- 

-1- 

+ + 


C5 

+ 

, , 


C6 


C7 


C8 

-1- 

-1- _^ 


Gl + 



0.0 0.2 0.4 0.6 0.8 1.00.0 0.2 0.4 0.6 0.8 1.00.0 0.2 0.4 0.6 0.8 1.( 

^debris/ ^tot ^debris/ ^tot ^debris/^tbt 



Figure 18. Separation of maxima and the fraction of debris par- 
ticles. Grey: 80% or more particles in debris, black: 20% or more 
particles are still bound. 



fitted a 'double Gaussian' to the energy distribution of tlie 
subhaloes in all our nine host haloes 



dN , - 
— = Aie 
aq 



+ A2e 2' 



(20) 



for a number of sample subhaloes are given by the solid black 
lines in Figure [161 

To gauge the integrity of our claim that it is possible to 
separate leading and trailing arm using the energy distribu- 
tion, we make use of the fact that we are in the unique sit- 
uation to separate leading and trailing arm particles again. 
We therefore calculate for each arm the energy distribution 
individually. We further decompose the fit to the combined 
energy distribution based upon all particles in the debris 
field (cf. Eq. (|20|l ) into two single Gaussian (i.e. one should 
describe the leading and the other the trailing arm). We 
check the validity of the decomposed Gaussians by compar- 
ing them to the corresponding "true" energy distributions 
derived via our "exact" leading/trailing arm classification 
method. Hence, for each subhalo we calculate two values 
representative of the squared deviations between the decom- 
posed fit and the true (binned) distributions for leading and 
trailing arm. The results can be viewed in Figures 1171 1181 



and 



191 to be explain ed in more detail below. 



where Ai,A2,Xi,X2,o-i, and CT2 are free parameters. The 
resulting best fit curves to the combined distribution dN/ dq 



Johnstoiil (| 19981 ) stated that the bimodality is in general 
more pronounced for subhaloes that are not disrupted, and 
for disrupted subhaloes a single-peaked distribution is ex- 
pected. We therefore start by inspecting the quality of our 
double-Gaussian fits by plotting the respective values 
against the fraction of mass in the debris field, i.e. number 
of particles ridebris in the debris field over number of particles 
ntot = 'T-dcbris + J^sat whcrc Tigat are the remnant particles in 
the subhalo; a ratio of ridcbris/ntot ~ 1 refers to a completely 
disrupted subhalo. The result is presented in Figure flTl where 
we observe a (marginal) trend for more disrupted subhaloes 
(i.e. larger ridebris/i-tot ratios) to hav e higher \^ valu es in 
agreement with our expectation and IjohnstonI (| 19981 ). re- 
spectively. However, we find no difference for leading (grey 
symbols) and trailing (black symbols) arm. 

In Figure [18] we plot the separation of the maxima of 
the two Gaussians (derived via fits to the distributions stem- 
ming from the separation method and not the decomposition 
of Eq. (|2U|I versus the debris fraction ndobris/Ti-tot again. The 
peak separation is further normalised to the sum of the half 
width at half maximum of each curve. We can confirm a 
small trend towards smaller separations for more disrupted 
subhaloes, i.e. the maxima are approaching each other and 
leading to less distinguishable curves for more disrupted sub- 
haloes, as expected. 

Finally we turn to our two subhalo properties again, 
namely the infall eccentricity and infall mass, examining if 
or how they relate to the (bi-modal) energy distribution. We 
plot in Figure [T5] the correlation of the infall eccentricity 
(left panel) and infall mass (right panel) with the quality of 
the decomposition of the fit of the energy distribution to a 
double Gaussian x^- We (again) observe little relation with 
eccentricity while there is a trend for more massive subhaloes 
to lead to better fits and hence a clearer separation between 
leading and trailing arm. Please note that we separated the 
subhaloes in Figure[T9]according to disrupted (grey symbols) 
and alive (black symbols) , where our criterion for disruption 
was fixed at ridebris /'T-tot ^0.8. 

6.6 Summary — Leading vs. Trailing Arms 

In this section we investigated trailing and leading arms in- 
dividually working out differences and similarities between 
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Figure 19. Correlation of tlie infall eccentricity (left panel) and infall mass (right panel) with the quality of the decomposition of the 
fit of the energy distribution into leading and trailing arm distribution. Grey (black) symbols mark disrupted (alive) subhaloes where 
the criterion for disruption corresponds to ndebris/"tot > 0.8 again. 



the debris fields. Somewhat surprisingly, most of the proper- 
ties we considered show little (if any) variation between the 
arms (e.g. velocity dispersion), but there does appear to be 
a mass difference between them. Our main conclusion is that 
even for complicated geometries of both arms it is in prin- 
ciple possible to separate individual particles/stars amongst 
leading or trailing arms by means of their energy distribu- 
tions. This in turn allows to (observationally) weigh both 
arms and verify our claim that their masses are different. 

In short, our results of the analysis of leading and trail- 
ing arm properties can be summarised as follows: 



• We presented a method to separate leading and trailing 
arm particles in cosmological simulations. 

• We confirmed that the method to separate leading and 
trailing arm utilising the energy distribution works well in 
cosmological simulations. 

• Trailing and leading arm can be separated observation- 
ally, when investigating the energy of debris particles: the 
energy distributions can be fitted by Gaussians, where the 
peak for the leading arm lies at smaller energies than the 
peak for the trailing arm. 

• The leading arm contains slightly more mass than the 
trailing arm for more than half of the subhaloes in our cluster 
haloes, whereas the opposite relation holds for the galactic 
halo. For more massive subhaloes, the fraction of mass in 
the trailing arm is increasingly higher. 

• We could not find differences in the velocity dispersions 
of leading and trailing arm. 



7 SUMMARY AND CONCLUSIONS 

In this paper, which is the first in a series, we have inves- 
tigated the tidal disruption of subhaloes (or satellite galax- 
ies) orbiting in cosmological dark matter host haloes. We 
used nine high resolution cosmological A'^-body simulations 
of individual galaxy- and cluster-mass haloes, each of which 
contain in excess of ~ 1 million particles within the virial 
radius, and within which a few hundred subhaloes can be 
reliably resolved. The time sampling of our simulations is 
such that we can track subhalo orbits in detail, which al- 
lows us to both determine orbital planes for subhaloes and 
to follow the mass loss they suffer as a function of time. 

In the first part of the paper (i.e. Section |3} we con- 
cluded that, while all of our host haloes are sufficiently re- 
laxed, they exhibit a range of integral properties. For exam- 
ple, they show a spread in formation times, mass assembly 
histories and shapes. These are precisely the kind of host 
properties that we expect to affect the properties of tidal 
debris field stripped from orbiting satellites. 

The second part of the paper (i.e. Sections |4]{6]) focused 
on the subhaloes and their tidal disruption. Our main con- 
clusions may be summarised as follows; 

• We predict that "backsplash" satellite galaxies should 
be found not only in the outskirts of galaxy clusters but also 
around galaxies. 

• Backsplash subhaloes contribute as much as 13% of the 
mass of the total tidal debris field within the host. This 
suggests that the region outside the virial radius of the host 
galaxy should be considered also when observing the debris 
of satellites and searching for possible corresponding satellite 
galaxies. 

• Tidal stream properties correlate well with progenitor 
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satellite properties. For example, we find a relationship be- 
tween the scatter about the (best-fit) debris plane and the 
infall mass of the satellite, as well as the spread in radial 
velocity and the infall eccentricity. 

• We devised a method to separate leading and trailing 
debris arm in cosmological simulations. 

• The masses of leading and trailing arms differ. 

• We do not find differences in the velocity dispersions of 
leading and trailing arms. 

• We present a method that uses energy distributions to 
separate leading and trailing arms and that can, in principle, 
be employed by observers. 

• The formation of the tidal debris field is extremely com- 
plex in "live" host haloes such that there appears to be little 
correlation between debris properties and host halo proper- 
ties (age, mass assembly history, dynamical state, shape). 

We set out to address two central questions key to un- 
derstanding tidal streams; 

• Can observations of tidal streams be used to infer the 
properties of their parent satellite? 

• Can tidal streams reveal properties of the underlying 
dark matter distribution of the host halo? 

These have proven to be very weighty questions to ad- 
dress, and so we focused on the first question in this paper. 
As the bullet points above indicate, the answer would appear 
to be yes. For example, there is a clear correlation between 
mass loss and orbital infall eccentricity, and between devia- 
tion of the tidal debris field from its orbital plane and the 
infall mass of the progenitor satellite. Addressing the second 
question is less straightforward. We have examined our data 
for correlations between stream properties and host prop- 
erties as summarised Table [T] but evidence for correlations 
with host mass accretion history, dynamical state and shape 
are extremely weak at best. Indeed, the only suggestion of a 
correlation that we find is between the mass ratio of leading 
to trailing arm and the virial mass of the host; while all our 
cluster mass haloes have more mass in the leading debris 
arm, the relation is reversed in the galaxy mass halo. 

Further, it is interesting to note that the galaxy halo 
Gl appears to be distinct from the eight galaxy cluster 
haloes (with the possible exception of CI) in terms of how 
its bulk properties relate to its stream properties. However, 
it is not clear whether this is a systematic effect, given 
that Gl is the only realisation of a galactic type halo. For 
example, the differences are not so striking if one compares 
Gl with CI. We are in the process of building a sample of 
galaxy haloes to supplement Gl. 

We now conclude with some general remarks. Previous 
studies of individual satellites disrupting in analytic host 
potentials using controlled A''-body experiments have hinted 
at simple and direct links between, for example, the flatten- 
ing of the host potential and the disper s ion of the debris 
field (e.g., Ilbata et al.ll2003l : lHelniill2004l : [Peiiarrubia et al] 
I2OO6I ). However, our experience suggests that the situation 
is vastly more complicated in "live" host haloes, where the 
host undergoes a complex mass accretion history, contains 
a wealth of subst ructure and cannot be easily modeled as 
a simple ellipsoid. iPefiarrubia et al] (|2006l ) pointed out that 
all present-day observable stream properties can only con- 



strain the present mass distribution of the host halo, in- 
dependent of its past evolution. We would argue strongly 
that even present-day properties of the host halo will be 
difficult to infer from stream properties, given the complex- 
ity of interactions responsible for the emergence of debris 
fields. There are multiple processes driving these effects and 
so while we note tentative correlations, we require further 
simulations to help us understand what is happening. This 
will represent the focus of the following paper in this series. 
We point out th at this situation is akin to the findings of 
I Gill et al] (|2004l ): they presented a detailed study of the dy- 
namics of satellite galaxies in similar host haloes and were 
also unable to establish correlations to the host properties. 

We close with a cautionary note: this paper is not to 
be understood as a definitive quantitative analysis of tidal 
streams in dark matter haloes. Rather it shows that cer- 
tain (unexpected) correlations between the debris field and 
satellite properties exist. Our study therefore represents the 
first step towards our understanding of whether or not tidal 
streams from disrupting satellite galaxies can be used to de- 
duce robustly properties of their host dark matter haloes 
and of the orbits of the satellites. 
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APPENDIX A: MEASURING HALO SHAPES: A 
NEW APPROACH 

Al Motivation for our Approach 

There is a rich literature on methods for the determination 
of halo shapes. Many of these studies utilise methods that 
identify subsets of particles and deduce axis ratios and ori- 
entations of best fitting ellipsoids to the haloes from the 
eigenvalues and eigenvector s of the moment of inertia ten- 
sor o f the particles (e.g., iFrenk et al.l 1 198^ : Ijing fc Sutd 
lioo^R However, the predicted shape depends on what 
is measure d. Typically parti cles are selected according to 
radius (e.g. iFrenk et al.lll988l ) b ut they can also be selected 
according t o local density (e.g. Ijing fc Sutd 1200 2| ) or po- 
tential (e.g. iPoweil I2OO3I ). with differing results. The pres- 
ence of sub structure within haloes al so affects the measured 
shape (e.g. iBailin fc Steinmet j[2005l '). We wish to develop a 
method that produces an unbiased estimate of a dark mat- 
ter halo's shape and orientation and that is unaffected by 
substructure. 

A2 Isodensity and Equipotential Surfaces 

It can be argued that the shape of the potential rather than 
the mass density provides a more meaningful measure of the 
shape of the halo for our purposes, because it is the potential 
that determines the orbits of satellites and tidal debris (e.g. 
iHavashi et al.1120071 ) . We use a novel technique to define both 
isodensity and equipotential surfaces while simultaneously 
allowing identifying a corresponding (spherical) radius Riso 
with the given surface. 

Our method for locating particles on isoden- 
sity /equipotential surfaces is as follows; 

At this po i nt w e need to men tion the novel met hod of 
ISpringel et~all i2004) (also used in iHavashi et"aIll2007^ which 
does not use the moment of inertia of the particle distribution to 
determine halo shape and orientation. Rather, it explicitly fits an 
ellipsoid to the potential estimated on a fine mesh. 
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• We start by calculating the local value of the density 
and potential at each particle position. To this extent we 
invoke again the original adaptive grid hierarchy used during 
the actual simulation, obtain the density and potential on 
these grids and interpolate to the particle positions. 

• We select a fiducial radius Riso- We then successively 
increase the thickness AT? of a spherical shell centred about 
the fiducial radius -Riso, i.e. the shell extends from -R™o" = 
Riso - AR to _R™ = Riso + AR. The value of AR is deter- 
mined in a way that we require to include a certain fraction 
of the total number of halo particles in the shell; in our case 
this amounts to 10%. 

• For all the particles in the shell we calculate the mean 
value of the density (potential) . This now defines the thresh- 
old isodensity (equipotential) value at our fiducial radius 

Riso- 

• From all simulation particles in the vicinity of the host 
(i.e. out to 2.5 i?vir) we now pick those ones whose local 
density (potential) value lies within a range of 10% about 
the mean value. 



A3 Shape Determination 

Having identified the subset of particles with the desired 
local density (potential) at 7?iao (cf. Section [A2|) . we com- 
pute their moment of inertia tensor Ijk- The eigenvalues and 
eigenvectors of Ijk are related to the axis ratios and orienta- 
tion, respectively, of the best-fitting ellipsoid to the particles. 
For a distribution of discrete point masses, Ijk is expressed 
as 

N 

Ijk ='^^mi{r'i&jk - XijXik) with j, fc = {1; 2; 3} , (Al) 

i = l 

where rrii is the mass of par ticle i, N the number of particles 
and ri = is the distance of particle i from 

the centre of mass of the particlef^. 

The novelty in our method of determining the shape is 
that we weigh particles by their local density when comput- 
ing Ijk as follows; 

JV ^ 

= ^■Z_. — {riSjk -XijXik) with j, k = {1; 2; 3}. (A3) 

i—i 

Here pi is the local density of particle i and the normalisation 
factor A is given by, 

1 



N 



(A4) 



This is necessary to ensure that especially particles in 
equipotential shells, which can have a substantial spread 
in their local densities, have approximately equal weighting 



in the summation. If this weighting is neglected, particles 
in higher density regions will dominate the summation and 
will bias our estimate of the axis ratios of the best-fitting 
ellipsoid. The formula for obtaining the correct axis lengths 
a > b > c from the eigenvalues of this inertia tensor is given 
in Appendix iBl below. 

As our host haloes contain a considerable amount of 
substructure which may distort the determination of the 
host's shape, we performed the shape analysis with as well 
as without removing all subhalos (i.e. all bound spherical 
halos found by the AHF analysis) from the host. We found 
that these subsystems do in fact have a great infiuence on 
the shape (especially on the shape of the potential) and thus 
we removed them for a credible and unbiased shape deter- 
mination. 

Our shape determination of all nine host haloes is pre- 
sented in Appendix[C]where we show the best fit ellipsoid to 
both the equipotential (upper panels) and isodensity (lower 
panels) contours at Riso ~ Rvir- 

We further calculate the usual measure for triaxiality 
T = (a^ - b^)/{a^ - c^) and list the values for both T and 
c/a (sphericity measure) at Tiiao = Rvii in Table [T] 



A3.1 Comparison to Previous Methods 

The simplest method determines the shape and orientation 
of a halo using all partic les within a spheri c al volume or shell 
at a given radius (e.g. Frenk et al. 19881: Kasun fc EvrardI 



at a given radms (e.g. l-brenk et al.l lineal: IKasun &: -bvrardi 
I2OO5I : iHopkins et al.|[2005l : iBailin fc Steinmetj|2005h . While 
this method robustly recovers the orientation of the halo, 
the resulting axis ratios tend to be biased towards larger 
values (i.e. haloes are predicted to be rounder). 

An alternative iterative approach to the problem again 
determines the shape and orientation of a halo using all 
particles within a spherical volume or shell, but this ini- 
tial surface or shell is now deformed in the direction 
along the principal axes of the best fitting ellipsoid; this 



Katdliggil: iDubinski fc Carlbers 


I1991I: IWarren et al.lll992l: 


BuUockl 200a Alkood et al.ll200e 


: Maccio et al.l 


2006b. Both 


Jins & Sutd ( 2OO2I) and iBailin & Steinmetd ( 


20051) have 



noted that iterative methods have difficulty in achieving 
convergence in simulations in which haloes are very well re- 
solved and contain a population of satellites. Satellites tend 
to lead to a distortion of the shape, and this is most pro- 
nounced in the outermost parts of host haloes where recently 
accreted satellites are most likely to be found. 

The impact of substructure can be reduced by working 
with a "reduced inertia tensor" ; 



M,k = 



E 



(A5) 



We note that the same ellipsoid fit can be obtained from the 
tensor 



M 



3 k 



N 

^ ^ m,iX{jXik 
i = l 



(A2) 



which has been used widely in previous studies. This differs from 
Ijk, which has an additional r?5jk- Both forms provide axis ratios 
and orientations that are identical (though the individual eigen- 
values are different). 



Here each particle is weighted by the inverse square of its dis- 
tance to the centre of the halo. While this recovers accurately 
the orientation of the ellipsoid, the axis ratios are systemati- 
cally overestimated and thus haloes are predicted to be more 
spher ical than they actually are (e.g. IBailin fc Steinmetd 
|2004 ). 

Both of these approaches have in general been applied 
in studies in which the best resolved haloes contain sev- 
eral tends of thousands of particles at most, in which it is 
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not possible to defined thin shells containing many particles 
with similar local densities or p otentials. A couple of studies 
(jjing fc Sutdl2002l : |Poweill2003h have looked at the shapes of 
well resolved haloes, each containing several hundred thou- 
sand particles within the virial radius. They identified shells 
of particles selected by local density and potential and deter- 
mined the best fitting ellipsoids. While this approach is sim- 
ilar to the one we have adopted, we note that these authors 
selected particles according to a local density or potential 
threshold. 



A3. 2 Advantages of Our New Approach 

The method presented in Appendix refappdsoequi enables 
us to select particles explicitly according to both their radius 
with respect to the centre of the halo and their local density 
(potential). This is desirable for a number of reasons. 

• It overcomes the problems of convergence faced by iter- 
ative methods and it avoids the artificial bias towards more 
spherical shapes introduced by the reduced inertia tensor 
feauation lA5|) . 

• By focusing on thin shells, we can obtain uncorrelated 
results for the shapes in the inner and outer parts of haloes. 

• By fixing the radius at which we select particles rather 
than fixin g the local density or potential threshold of parti- 
cles (as in ljing fc Sutoll2002l : |Poweill2003D . we can compare 
directly the shape of the local density and potential at the 
same radius. 

• We can compare the shapes of different haloes at the 
same fraction of the virial radius, which depends solely on 
virial mass, rather than at a similar local density or po- 
tential, which will depend on concentration, recent merging 
history, larger scale environment, etc.... 



APPENDIX B: MEASURING HALO SHAPES: 
FROM EIGENVALUES TO ELLIPSOID AXES 

When we diagonalise the moment of inertia tensor (as de- 
fined by equation [XT} for a thin shell of particles, we obtain 
the "principal moments". 



i = l 
N 

A2 = ly'y' = My^(2:f + X'^) 



1 = 1 

AT 



A3 = i-j,j^MY,{^? + y?)- 



(Bl) 
(B2) 
(B3) 



Here x ,y , z are the coordinates in the principal axes frame, 
M is the total mass of the shell. 

For a thin ellipsoidal shell with major, intermediate and 
minor axes a > b > c, the corresponding principal moments 
are. 



Acii.i = ^M{b^ + c^) 
Aoii,2 = ^M{c^+a^) 
Aoii,3 = ^M{a^ + b'') 



(B4) 
(B5) 
(B6) 



If the thin shell of particles can be described as an ellipsoidal 
shell, then the coordinate axes a, b and c can be chosen such 
that Acii,i = Ai. It follows that 



3(-Ai -h A2 A3) 



2M 



b = 



3(Ai - A2 + A3) 



2M 



3(Ai + A2 - A3) 



2M 



(B7) 
(B8) 
(B9) 



which allow us to relate our thin shells of particles to their 
best fitting ellipsoids. 

This paper has been typeset from a file prepared 

by the author. 



APPENDIX C: HOST HALO SHAPES 
(C1-C8 AND Gl) 

Here we present our shape determination of all nine host 
haloes by showing the best fit ellipsoid to both the equipo- 
tential (upper panels) and isodensity (lower panels) contours 

at -Riso ~ -Rvir . 

Even though most ellipsoids fit rather well, we also note 
that for some host haloes there is an offset with respect to 
the particle distribution (especially C6, C8, Gl). In these 
cases, the density centre of the host does not correspond to 
the centre of isodensity (equipotential) shells. In most these 
cases (i.e. C6 and Gl) there is a more or less pronounced 
'bump' visible in the equipotential plots (upper panels), indi- 
cating an adjacent substructure currently merging with the 
host halo. Though substructure was removed, small over- 
densities in the environment of the subhalo still r e main and 
distort the equipotential surfaces. iGao fc White! (|2006l ) re- 
cently studied asymmetries in the inner regions of ACDM 
haloes, claiming that recent accretion events and deviations 
from equilibrium are responsible for such an offset. 

We further note that the isodensity as well as the 
equipotential surfaces in Figures rCll - IC9l are in fact defining 
halo boundaries close to the virial radius (dashed circles), 
just as intended. Or in other words, our method ensures 
measuring the shapes of the haloes close to the pre-selected 
_Riso value. 



APPENDIX D: FINDING THE APPROPRIATE 
TUBE RADIUS AROUND THE ORBIT 

Dl Iteration 

For a given tube radius, we count all particles lying inside 
the tube around the orbit, thus determining the fraction of 
'tube particles'. If more than 68.3% (1 a) of the particles 
are lying inside, the radius is decreased, otherwise enlarged. 
The procedure is repeated, until we achieve a tube particle 
fraction of 68.3% ± 0.5%. If no convergence is reached, we 
increase the tolerance by a factor of 2, up to 2%. In cases of 
deviations exceeding 2%, we skip the satellites. This usually 
is the case, if the number of debris particles is so small, that 
a meaningful statistic is not possible anymore. 
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Figure CI. Best fit eUipsoids (solid lines) for halo CI at 
Riso=RviT- The top row shows particles in the equipotential shell, 
projected to the planes of the ellipsoid (eigenvector coordinate 
system). The bottom row shows the best fit ellipsoid for the cor- 
responding isodensity range. The dashed circles mark the virial 
radius. Particles bound to satellites have been removed and re- 
maining particles weighted by their local density prior to fitting. 



Figure C4. Ellipsoid fit for the cluster C4. 
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Figure C2. Ellipsoid fit for the cluster C2. 
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Figure C5. Ellipsoid fit for the cluster C5. 




Figure C6. Ellipsoid fit for the cluster C6. 



Figure C3. Ellipsoid fit for the cluster C3. 
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radius rt around the orbit path (thick black line). Not only a 
cylinder around the orbit path is used, but also the region in a 
Figure C7. Ellipsoid fit for the cluster C7. half sphere at the end of each orbit piece is included. 




Figure C8. Ellipsoid fit for the cluster C8. 
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D2 Particles inside/outside the tube 

We check for each particle, whether it Ues inside a tube of a 
given radius or outside. 

Instead of first searching for the orbit point closest to 
the considered particle, we go through the orbit points until 
one is found, where the particle lies in the tube next to it. 
Thus, we find one region where the satellite particle lies in- 
side the tube. In fact, there could be more than one suitable 
region, maybe at some other point the particle would even 
lie closer to the orbit path. Yet knowing that the particle lies 
inside the tube at some point is sufficient for our purposes. 

Figure IDll illustrates the situation of particle P lying 
outside the orbit tube. A particle at distance Vp from an 
orbit point lies inside the corresponding part of the orbit 
tube with radius rtubc, if 

• '"p,± < '"tube (within cylinder) 

• '"p.ii ^ ftubo (include cap of radius rtubo behind the 
considered orbit point) 

• '"p,|| ^ \ ^ n| +'"tubo (not too far beyond next orbit 
point) 

with the orbit points n and n+i and the components of Vp 
parallel (rpji) and perpendicular (fp^x) to the orbit path. 
Particles are not only checked for lying in the cylinder of 
radius rtubc around the orbit path, but also a small half 
sphere ('cap') is allowed on front of and behind every orbit 
point. Otherwise, particles at the kink of the orbit might be 
wrongly classified as non-tube particles. 



Figure C9. Ellipsoid fit for the galaxy Gl. 



